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' Abstract. We describe similariton pulse propagation in double-doped optical fibers with the 

CN ' aid of self-similarity analysis of the cubic-quintic nonlinear Schrodinger equation with varying 

dispersion, nonhnearity, gain or absorption, and nonhnear gain. Exact similariton pulses that 
can propagate self similarly subject to simple scaling rules of this model have been found using 
a fractional transform. By appropriately tailoring the dispersion profile and nonhnearity, the 
condition for optimal pulse compression has been obtained. Also, nonlinear chirping of the 
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1. Introduction 

Over three decades ago [1], a major contribution to fiber-optics research was made when 
Hasegawa and Tappert pointed out that soHtons could propagate in glass fibers with anomalous 
group velocity dispersion (GVD) under the influence of the optical Kerr effect. Optical 
solitons were observed in 1980 [2], and are of fundamental interest because of advances 
in the control and generation of ultrashort optical pulses [3] and their potential applications 
in the field of optical fiber communications. Recently, solitary wave pulses having linear 
chirp and parabolic intensity profile have been predicted and experimentally observed in 
an optical fiber [4] and in optical fiber amplifiers [5, 6] with normal GVD. These optical 
pulses propagating in a self-similar regime were called similariton pulses, or more pertinently 
parabolic similariton pulses (PSPs). These pulses have recently found increasing practical 
application in high-power amplifier systems [7, 8, 9, 10], efficient temporal compressors 
[5, 7, 8], and similariton lasers [11]. Another type of solitary wave pulse propagating in 
self-similar regimes with linear chirp and sech-amplitude profile has been predicted in optical 
fiber amplifier with anomalous GVD [12, 13, 14, 15, 16]. These pulses, for the same reason 
were called similariton pulses, and to distinguish them from PSPs we will also call them sech- 
similariton pulses (SSPs). These SSPs propagating in self-similar regimes are of considerable 
interest for pulse compression applications that recently were demonstrated experimentally 
[17]. 

By now, it has been well established that cubic-quintic nonlinear Schrodinger equation 
(CQNLSE) with nonlinearity management presents practical interest since it appears in 
many branches of physics such as nonlinear optics and Bose-Einstein condensate (BEC). 
In nonlinear optics it describes the propagation of pulses in double-doped optical fibers 
[18]. In BEC it models the condensate with two and three body interactions [19, 20]. In 
optical fibers periodic variation of the nonlinearity can be achieved by varying the type 
of dopants along the fiber In BEC the variation of the atomic scattering length by the 
Feshbach resonance technique leads to the oscillations of the mean field cubic nonlinearity 
[21]. Serkin et al described optical organic materials and waveguides modeled by CQNLSE 
as the key elements for future telecommunications and photonic technologies [22]. For 
example thin films of polydiacetylene para toluene sulfonate exhibit cubic-quintic nonlinearity 
[22]. Modulational instability of the CQNLSE through the variational approach has been 
discussed [23]. Recently, solitons with cubic and quintic nonlinearities modulated in space 
and time has been discussed [24]. An analytical approach to soliton of the saturable nonlinear 
Schrodinger equation determination and consideration of stability of solitary solutions of 
CQNLSE has been discussed [25]. Also, solitary wave solutions for high dispersive CQNLSE 
has been discussed [26]. Although we discuss self-similar autonomous solutions in this 
model, interested reader is referred to excellent reviews by Serkin et al in Refs.[27, 28] on 
nonautonomous solitons. In this paper, we describe a different chirped similariton pulses 
for the generalized CQNLSE with varying dispersion, nonlinearity, gain or absorption, 
and nonlinear gain using a recently developed fractional transform (FT) [29]. Unlike the 
PSPs and SSPs, the similariton pulses found for this model are Lorentzian-type ones some 
being singular, and others being nonsingular, and periodic. Furthermore, we delineate the 
compression of these exact chirped similariton pulses by piloting, precisely, the pulse position 
through appropriate tailoring of the dispersion profile. The results reported here have potential 
application to the design of double-doped fiber optic amplifiers, optical pulse compressors, 
and solitary wave based communication links using higher order nonlinearity and nonlinear 
gain. 
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2. SIMILARITON PROPAGATION IN DOUBLE-DOPED FIBERS 

As detailed in Ref. [18], the double-doped semiconductor fiber is fabricated on the basis of 
silica glass doped by two appropriate semiconductor compounds. In that case the nonlinear 
correction to the medium's refractive index can be expressed in the form 6n = nxl - n^I^, 
where / being the light intensity and the coefficients n2 > and n4 > determine the 
nonlinear response of the media. They are related to third order susceptibility x^^^ and fifth 
order susceptibility ;^''^^' through n2 - ?>x^^^ l^na and n^ - -5x^^^ l^'^no, where no being the 
linear refractive index. Although, formally it may be obtained by an expansion of the saturable 
nonlinearity 6n - n2l{\ + (n4/n2)/]"', it is restricted under the effect of self- focusing as 
d{6n) - dl. 

The generalized CQNLSE with distributed nonlinear gain governing the propagation of 
the optical field in a single-mode double-doped optical fiber can be written in the form 

i<Pz = ^tArr - r(z) I lA I V + 6(z) I <A iV + '9(z)'A + ixiz) I "A I V, (1) 

where if/ is the complex envelope of the electric field in a comoving frame, z is the propagation 
distance, t is the retarded time, g{z) is the gain function, and ;^f(^) accounts for the nonlinear 
gain or absorption [30, 31]. In the absence of the last term and quintic nonlinear term, this 
equation has exact similariton solutions that exhibit linear chirp. So also, in the presence 
of nonlinear gain and in the absence of quintic nonlinear term, this equation has exact 
chirped SSPs. Recently, pedestal free compression of the pulses in the generalized nonlinear 
Schrodinger equation with quintic nonlinearity without nonlinear gain term has been reported 
[32]. Also, chirped and chirp free self-similar cnoidal and solitary wave solutions of this 
model without nonlinear gain term have been analyzed [33]. But here we are concerned with 
similariton pulses characterized by a nonlinear chirp, resulting from the nonlinear gain. For 
finding the similariton solutions using the FT of Eq. (1), one writes the complex function 
il/{z, t) as 

iA(z, t) = P{z, T)exp{/moln[P(z, t)] + i^(z, t)) (2) 

where mo denotes the nonlinear chirp parameter, and P and O are real functions of z and t. 
Using this ansatz, we find the system of two equations for phase <l)(z, t), and amplitude P(z, t) 
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In general case when the coefficients of the CQNLSE with gain and nonlinear gain are the 
functions of the distance z, the amplitude P{z, t) of the similariton solutions has the form 

Q(T) 
P(z, r) = -^exp(5 (z) - mo0(z)), (5) 

where the scaling variable T is given by 

T = ^^^^. (6) 

r(z) 

And the other functions Riz), &iz), S (z) and Tp(z) take the forms 

R(z)^ Cf^izydz', (7) 
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where bo, A, cq and t^. are the integration constants and p(z) = /3(z)/y(z)- Here T(z) and Q{T) 
are some functions we seek, where without loss of generahty we can assume r(0) - 1 . Let us 
assume that the phase has the quadratic form 

<5(z, r) = a(z) + b(z)[T - Tpiz)] + c(z)(t - t„)\ (11) 

where the pulse position t^ is a function of z. Thus we consider the class of similariton 
solutions with the phase given by Eq. (1 1). Then Eq. (3) with the phase (1 1) can be written 
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where 



da Idb dT„(z)\ drJz) dc(z) -, 

M = — + — - 2c(z)^— (T - Tpiz)) - b(z)^— + -f^(T - Tp(z)fXl3) 
dz \dz dz dz dz 



N^P [b(z) + 2C(Z)(T - Tp(z))f + 2moPr[b(z) + 2C(Z)(T - Tpiz))] + ^ - Prr- (14) 

Equation (12) contains an explicit dependence on the variable (r-Tp) which disappears when 
the terms at monomial (t - t^)^ are equals, hence we find three equations 
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Combining Eqs. (4) and (11) we find a third equation for our general class of similariton 
autonomous solutions with quadratic phase 

Biz) \ P^ 

P, = l3{z)c{z)P +/3{z)b(z)Pr + 2/3{z)c{z)(t - t^(z))P, + ^ moP,, + mo-^ 

Taking into account Eq. (5) we find that Eq. (19) will be satisfied if and only if the function 
r(z) is defined as 

1 dr(z) 



nz) dz --'^''^'''^- 

The solutions of Eqs. (15), (16), and (20) are 
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bo dR boP 

Pr = — P - -^, (22) 

mo/3 dz mo 

where b(z) = bo, and 

r(z) = 1 - 2cqR(z). (23) 

Taking into account Eqs. (5) and (17) we find 

d^Q 2yT , 26 , 

-l^-yQ- ^exp[25 - 2mo0]e' + — exp[45 - Amo(d]Q^ = 0. (24) 



dT^ "^ 13 "^^ '"^ p 



where 
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V = T^—^{dRldzf - r^—^(dR/dz) + F^— - T^bl - T^-(daldz) - 
mo/3^ mo/3 wo /3 

moT dr 2mor^ dS d& 

/3 dz /3 dz dz 

In general case, the coefficients in Eq. (24) are functions of variable z but the function Q{T) 

depends only on the scaling variable T, hence this equation has nontrivial solutions [2(7") ^ 0] 

if and only if the coefficients in Eq. (24) are constants 

b^ 2b^ b^ 2 

A = T^—^(dR/dzf - T^—^(dR/dz) + T^— - T^bl - r^-(da/dz) - 
mo/3'^ mo/3 mo /3 

moT dV 2moT^ dS d@ 

/3 dz /3 dz dz ' 

yr 

— exp[25 - 2mo&] - a, (27) 

^exp[45 - 4mo©] = e. (28) 

Here A - const, a = const, hence Eqs. (26) and (27) yield 

,2, 1-^0 , 2 da ml -2 

mo p(0) dz 2 

because r(0) = 1 and 5(0) = ^l"! ~l~P(0) I. Thus in nontrivial case Eq. (24) can be written 
as 

^-AQ-2aQ'+2eQ'^0 (30) 

The solution of Eq. (26) is 



I +ml bl mn 

a(z) = flo + —^®(z) - -fR(z) - Y^n[\ c(z) |] - moS(z) (31) 

It is useful to write Eq. (27) in the form 

p(z) = p(0)[l - 2co/?(z)]exp[25 - 2mo0] (32) 

where we define the function p(z) as 

,, m ... AO) 1 ,„, 

p(z) = — -; p(0) = -77; = -■ (33) 

riz) r(0) a 

One may differentiate this equation and find 

.. 1 d co/3(z) nAmo/3(z) 

'^'^ = MTUz'^'^ " T^^^^ - il-2coRiz)r ^''^ 
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^ = 4^ (35) 

y{z) ml -2 

where /i = 1 or |. The former condition describes that the parameter functions in Eq. (1) 
can not be chosen independently, while the latter condition implies that the nonlinear chirp 
parameter mo in fact be determined by the ratio xidlyiz). In this paper, it requires that this 
ratio is a constant. From the physical cases we have considered, we come to the conclusion 
that m^ ^ 2 for arbitrary nonlinear materials. 

3. Exact similar iton solutions 

In order to obtain the exact similariton solutions of Eq. (1), we use a fractional transform 

^ ' \+Dp(T) ^ ' 

that connects Eq. (30) to the elliptic equation: /" + af + bf^ - 0, where a and b are real. 
As is well known, / can be taken as any of the three Jacobian elliptic functions with an 
appropriate modulus parameter: cn(r, m), dn(r, m), and sn(r, m), with amplitude and width, 
appropriately depending on m. Using the limiting conditions of the modulus parameter, one 
can obtain both localized and trigonometric solutions. We list below a few localized as well 
as periodic solitary wave solutions. 

For explicitness, we consider Eq. (30), with all the parameters and illustrate below 
various types of solutions of Eq. (1), by taking / - cn{T,m). Other cases can be 
deduced by using the Jacobian elliptic function identities: dn (T, m) - 1 - msrp-(T, m) and 
sn^(r, m) - \ - cv^ij, m). The consistency conditions are given by, 

-AA + 2(AD - B)(m - 1) - 2aA^ + 2eA^ = 0, (37) 

2(AD - B)[2(l - 2m) + 3D(1 - m)] + 4-D(AD - B)(m - 1) - A(2AD + B)- 

2a(3A^B + 2A^D) + 10eA'*B = 0, (38) 

2(AD - B)[2D(2m - 1) + 3m] + 4-D(AD - B)[2(l - 2m) + 3D(1 - m)] + 2D\AD - B)(m - 1) - 

A(6AD^ + 2BD) - 2a{'iAB^ + 6A^BD + A^D^) + 20eA^B^ = 0, (39) 

-2mD(AD -B) + 4-D(AD - B)[2D(2m - 1) + 3m] + 2D^(AD - B)[2{\ - 2m) + 3D(1 - m)] - 

A{AAD^ + 6BD^) - 2a(B^ + 6AB^D + 3A^BD^) + 20€A^B^ = 0, (40) 

-4mD^(AD -B) + 2D^{AD - B)[2D{2m - 1) + 3m] - A{AD'^ + BD^) - 

2a(2B^D + 3AB^D^) + lOeAB* = 0, (41) 

-2mD\AD -B)- ABD'' - 2aB^D^ + 2eB^ = 0. (42) 

The above equations clearly indicate that the solutions, for m - 1, m = and other values of 
m, have distinct properties. To find the solution for either m - 1 or m = or other values of 
m we should keep all the parameters in the FT. Here, we list three such cases: m - 0,m - 1, 
and m - 1 /2 in detail. 

3.1. Trigonometric solution 

For m = 0, from Eq. (42) we find that 

B = VrZ) (43) 

where 

_ a+ Va^ + 2eA 
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From Eq. (41) we find that 

A = -!-^ -. — (44) 

A - (IQer^ - 6aY' - A) 

And we find that 

_ (4 + 2/1) VF - (-2a + 20er')A^ + \2aA^ y/T - (-6A - GaV + 4)A 

LJ y-TU ) 

ioa-ioVf 

Hence, the trigonometric solution of Eq. (1) can be written as 



,, , M-2mz)lyiz))A-.Bco.-^[l^^^] 

Uriz t) — 

^/2[l-2coR(z)] 1 + Dcos^ [l^I^i^j 



X expi/moln 



1 +Dcos2[l^i^i^2^1 



+ i<t>(z,T)}. (46) 



We emphasize that this trigonometric solution is the general solution of this model and is valid 
for all values of the amplitude parameters A, B, and D subject to the condition AD - B ^ 0. 
This solution of Eq. (1) has no analogue in the limit n4 — > 0. 

3.2. Hyperbolic solution 

For m-\, from Eq. (37) we are able to determine the amplitude parameter A completely as 

9 a + Va^ + 2eA 

A^ = . (47) 

2e 

From Eq. (38) we find that B and D are related to each other as 

B = FzD, (48) 

where 

A{AaA^ + 2A + A) 

r2 - 



lOeA^ - 6aA^ - A + A' 
Using this value of B in Eqs. (39) and (40) we fully determine D as 

n- ^iz^I^ ^. 

4r2[i + a(A + r2)]-8A-r2-^ 

Thus, the solitary wave solution of Eq. (1) is written as 

^(m2-2)08(z)/r(z))A + Bsech^ [ ^I^.^zf ] 



^{z,t)^ 



X exp I /mo In 



V2 [1 - 2cMz)] 1 + Dsech^ [ 7:f,Xf ] 



+ i^(z,T)\. (50) 



Again, we emphasize that this solution of Eq. (1) has no analogue in the limit «4 — > 0. 
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Figure 1. Density plot depicting the nonlinear chirping of the trigonometric solution for 
o- = 0.5, a = 3.45, A = -3.34, and e = -2.78 



3.3. Pure cnoidal solution 

In order to obtain pure cnoidal wave solution, we put m =1/2 and B - Din the consistency 
conditions. From Eq. (40) we determine the amplitude parameter A completely as 



-b± V/?^ - Aac 
2a 



(51) 



where a - 20e - 6a, b -3 -4A- 12q;, and c - 10(e - a) - llA-2a -3. We determine the 
value of D by solving Eq. (42) as 

A- 1 



2(e - a)- A 
Hence, the pure cnoidal wave solution of Eq. (1) with all the amplitude parameters is given 
by 



<P(z,t) 



^{ml-2)(Ji{z)ly{z))A + Ben"" [ 



T-T,+boR{z) 
\-2caR(z) ' 



'] 



\lu 



V2 [1 - 2coR{z)] 1 + Z)cn2 [^yS^^ '«] 



X exp I HWoln 



^^b-A't^^M 



+ /<[)(z, t) 



.l+^cn2[l,I^t^,m] 
Once again, we emphasize that this solution of Eq. (1) has no analogue in the limit n^ — > 



(53) 



4. Nonlinear chirping 

In this subsection, we wish to cite an example corresponding to the trigonometric solution, 
illustrative of the fascinating features of chirping in double-doped fibers, by considering the 
system in which the GVD and the nonlinearity are distributed according to 

j8(z) = Po cos(crz), y(z) = Jq cos(crz), (54) 

where ySo, To and a-{i^ 0) are arbitrary constants. In this case, the corresponding gain, and the 
nonlinear gain of the fiber amplifier are given by 

crv cos(crz) yU/lmo cos(crz) 



0(Z) = 



2 - 2 V sin(cr2) ( 1 - v sin((Tz))2 ' 



(55) 
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Figure 2. Density plot depicting the nonlinear compression of the trigonometric solution for 
o- = 1.0, a = 3.45, A = -3.34, and e = -2.78 
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cos(crz). 



(56) 



where the parameter v = 2co/3o/o- has been introduced for brevity. 
Hence the ampHtude of the similariton pulse is 



P(Z,T) 



where 



V2(l - vsino-z) 1 + Dcos2 [^^j 



W(z) = To[l - vsin(crz)], 



(57) 



and the pulse position Tp varies with Tp = Tc-(bo/3o/cr)sm{crz)- The resultant chirp consisting 
of linear and nonlinear contributions are derived as [5] 



6u(t) 



mo 



tan(9|) 



W{z)l+DcosHg^) 



■ W(z) > 



2coTo 

W(z) ^' 



(58) 



We notice that the first term in Eq. (58) denotes the nonUnear chirp that results from the 
nonlinear gain, while the last two terms account for the linear chirp. The propagation of this 
chirped pulse has been depicted in Fig. 1 for various parameter values of tq and v = 1 . 



5. Nonlinear compression 

Next, we elucidate the compression problem of the pulse in a dispersion decreasing optical 
fiber. For the purpose of comparison with Ref. [14, 15], we assume that the GVD and the 
nonlinearity are distributed according to the following relations 

j0(z) = /3o exp(-o-z), y(z) = yo exp(Q'z), (59) 

where j6o < 0, yo ^ 0, and cr + ^. Then from Eq. (34) the expression for gain can be 
calculated as 



giz) 



-a + vcr- 



(60) 



1 - v( 1 - e-'^-) ( 1 - v( 1 - e-'^'))2 ■ 
Let us consider the most typical physical situation when the loss in an optical fiber is 
constant i.e., y{z) is constant. According to Eq. (60), this occurs when v = 1 and a > 0; 
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g(z) = -a, hence the gain is negative. It is remarkable that the width of the solutions presented 
here tend to zero when z — > oo. 

We apply the above insights to nonlinear compression of the trigonometric solution given 
by Eq. (46). Then the amplitude of the compressed wave is 

A + Bcos2 [I^f 1 

1 + Dcos2 [i^l 

where 



J(ml-2)p(0)exp[-(cT + a)z] 

Uiz) = ^^ . 

V2[l-v(l-exp(-crz))] 

And W(z) = To[l - v(l - exp(-(Tz))]. Fig. 2 shows that for constant loss the trigonometric 
pulse can be compressed to any required degree as z — » oo, while maintaining its respective 
original shape. 

6. Conclusion 

In conclusion, for generalized cubic-quintic nonlinear Schrodinger-type equation with 
variable dispersion, variable Kerr nonlinearity, variable gain or loss, and nonlinear gain, we 
found exact chirped pulses that can propagate self-similarly subject to simple scaling rules, 
of this model. The fact that the pulse position of these chirped pulses can be precisely 
piloted by appropriately tailoring the dispersion profile, is profitably exploited to achieve 
optimal pulse compression of these newly reported chirped similariton solutions. Studying 
soliton bistability using the exact solutions found here will be an interesting issue to be 
pursued. These analytical findings suggest potential applications in areas such as optical fiber 
compressors, optical fiber amplifiers, nonlinear optical switches, and optical communications. 
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